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ABSTRACT 

This  report  deals  with  design  optimization  of  Internal- 
Expanding  Rim  Brakes.  A  computer  program  was  developed  to 
calculate  the  actuating  force,  torque,  stopping  time  and 
drum  temperature.  The  drum  temperature  is  calculated  by  the 
finite  difference  method. 

A  comparison  of  results  has  been  made  using  a  simplified 
equation  that  is  in  common  use  in  engineering  texts. 

Numerical  optimization  is  shown  to  be  a  convenient  tool 
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a  Distance   from    pivot    to   the    center   of   rotation    (m) . 
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b  Width    of   friction   laterial    (in)  . 

B.  Biot    modulus. 

1 

c  Specific  heat  (J/Kg-°C). 

C  Thermal  capacity  (J/°C). 

d  Distance  from  actuating  force  to  the  hinged  pin  (m) . 

dc  Rats  of  deceleration  (a/sec?  ) . 

E  Kinetic  energy  (J)  . 

f  Frictional  force  (N) 

F  Actuating  force  (N)  . 

F0  Fourier  modulus. 

g  Gravity  constant  (m/sec2.  )  . 

h  Convection  heat  transfer  coefficient  (W/m2-C). 
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M  Normal  moment  (N-m)  . 
n 
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th  \      f     t 
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T  Temperature  (°C)  . 

T0  Torque  (N-m)  . 

V  Velocity  (m/sec). 

V0  Volume  (m3)  . 

W  Vehicle  weight  (N)  . 


NOTAriON 

R. .    The  thermal   resistance    between    ncde   i    and   the 
iJ 
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I.   INTRODUCTION 

Brakes  are  mechanical  devices  for  retarding  the  motion 
of  a  vehicle  or  machine  by  means  of  friction.  Because  of  the 
similarity  of  their  functions,  many  clutches  may  also  be 
included  here,  assuming  centrifugal  forces  are  accounted 
for. 

A  simplified  dynamic  representation  of  a  brake  is  shown 
in  Fig.  1.  Two  masses  with  inertias,  I:  and  I2  ,  rotating 
at  the  respective  angular  velocities  co i  and  w2  (one  of  which 
may  be  zero) ,  are  to  be  brought  to  the  same  speed  by 
engaging  the  brake. 

The  friction  arake  has  three  basic  elements;  two 
opposing  friction  surfaces  and  a  mechanism  for  forcing  the 
friction  surfaces  into  contact,  whenever  a  friction  brake  is 
engaged  to  join  two  members  having  relative  motion,  there  is 
a  period  of  slip  which  may  last  several  seconds.  This  slip 
is  one  of  the  chief  merits  of  the  friction  brake;  i~  absorbs 
siocks  and  prevents  excessive  torsional  stresses  on  the 
power  transmission  system.  On  the  other  hand,  slip  is  the 
limiting  factor  in  friction  clutch  and  brake  performance; 
for  heat  is  generated  in  proportion  to  slip,  torgue 
transmitted,  and  period  of  slip. 

The  following  parameters  ars  of  interest  in  analyzing 
the  performance  of  these  devices; 

1.  The  actuating  force. 

2.  The  torque  transmitted. 

3.  The  temperature  rise. 

4.  The  slip  time. 

This  report  deals  with  Int arnal-Expanding  Him  Brakes. 
This  formulation  also  applies  to  internal-expanding  clutches 
if  centrifugal  forces  are  accounted  for. 
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II. 


INTERNAL-EXPANDING    HIM    CLUTCHES    AND    BRAKES 


A.  GENERAL    MACHANICAL    PRINCIPALS 

A  brake  or  clutch  assembly,  uses  a  brake  shoe  to  which 
is  attached  a  friction  material,  called  lining.  The  lining 
is  riveted  or  bonded  to  the  brake  shoe  as  shown  in  Fig.  2. 
The  brake  shoe  is  pivoted  at  a  fixed  point  and  the  other  end 
is  subjected  to  a  force  which  presses  the  shoe  in  contact 
with  the  drum.  The  force  between  the  brake  and  the  drum  is 
radial  as  the  drum  rotates.  If  a  point  on  the  rotating  drum 
surface  first  makes  contact  with  the  shoe  at  the  end  nearest 
the  pivot,  the  shoe  is  termed  a  "trailing  shoe".  If  it  first 
makes  contact  at  the  other  end  the  shoe  is  termed  "leading 
shoe",  the  latter  giving  a  higher  braking  torque  than  the 
former  for  a  given  braking  force. 

The  friction  between  the  lining  and  the  drum  creates 
heat  which  is  basically  the  conversion  of  energy  of  motion 
of  the  vehicle  or  machine  to  thermal  energy  at  the  friction 
surfaces,  namely  the  lining  and  the  drum.  This  heat  is  then 
dissipated  and  absorbed  by  the  drum  by  conduction, 
convection  and  radiation  into  the  atmosphere. 

B.  FRICTION  FUNDAMENTALS  AND  MATERIALS 

Friction  mechanisms,  such  as  brakes,  are  systems  for 
converting  mechanical  energy  into  heat.  Several  basic 
factors  affect  friction  and  wear  of  materials  used  in  brake 
systems.  The  main  factors  are  temperature,  pressure,  speed, 
surface  roughness,  and  type  of  material.  Some  organic  or 
molded  friction  materials  show  no  change  in  friction 
characteristics  with  pressure,  while  others  such  as 
sintered-metal  materials  decrease  in  friction  coefficient  as 
pressure  is  increased.  For  metallic  friction  materials  there 
is   also  a   decrease   in  coefficient   of   friction  as   speed 
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increases.  Temperature  effects  upon  the  coefficient  of 
friction  vary  widely  with  the  type  of  materials  used. 

In  a  two-shoe  internal  expanding  brake  there  is  a 
tendency  for  the  brake  drum  to  deform  under  hard 
application.  Drums  become  elliptical  and  the  force  to  do 
this  is  guite  high  and  contributes  to  friction  force. 

A  brake  or  clutch  friction  material  should  have  the 
following  characteristics  to  a  degree  which  is  dependent 
upon  the  severity  of  the  service: 

1.  A  high  and  uniform  coefficient  of  friction. 

2.  The  ability  to  withstand  high  temperatures,  together 
with  good  heat  conductivity. 

3.  Properties  whioh  are  not  affected  by  environmental 
conditions  such  as  moisture. 

4.  Good  resiliency. 

5.  High  resistance  to  wear,  scoring  and  galling. 

C.   BRAKE  DRUMS 

One  of  the  primary  functions  of  a  brake  drum  is  that  of 
absorbing  and  dissipating  the  heat  developed  during  the 
application  of  the  brake.  A  brake  drum  is  a  heat  sink  into 
which  heat  goes  after  it  is  created  by  the  rubbing  friction 
of  the  brake  lining  contact  to  drum.  The  brake  shoe  and 
lining  permanently  fixed  on  the  axle,  when  actuated, 
contacts  the  drum  under  pressure  to  cause  the  friction  to 
stop  the  vehicle.  The  energy  of  motion  of  a  vehicle  is 
converted  to  thermal  energy  by  the  brake  assemblies.  A  brake 
drum  must  have  the  capacity  to  absorb  and  dissipate  this 
hsat  energy  within  the  limits  of  the  brake  heat  input.  If 
this  is  not  the  case,  the  drum  expands  and  the  brakes  fade 
or  fail.  The  greater  the  mass  of  the  drum,  the  more  heat  it 
can  absorb  and  store  until  such  time  as  the  heat  can  be 
dissipated  by  convection  and  radiation  [Ref.  1]. 

An  ideal  brake  drum  would  have  the  following 
characteristics; 
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1.  High  structural  strength  to  resist  bursting  forces. 

2.  Uniform  coefficient  of  friction. 

3.  Hard  surface  to  resist  scoring. 

4.  High  heat  conductivity  to  rapidly  conduct  heat  away 
from  braking  surfaces. 

5.  High  emissivity  factor  to  radiate  heat  from  the  drum 
surface  to  the  atmosphere. 

6.  High  heat  storage  capacity  to  store  heat  from 
successive  brake  applications  until  it  can  be 
dissipated. 

7.  Good  nachinability  to  permit  boring  of  the  drum. 

D.   STATIC  AND  DYNAMIC  ANALYSIS 

1 .  Assumptions 

In  developing  the  equations,  the  following 
assumptions  have  been  made; 

a.  The  pressure  at  any  point  on  the  shoe  is  proportional 
to  the  moment  arm  of  this  point  from  the  pivot. 

b.  The  effect  of  centrifugal  force  may  be  neglected. 

c.  The  shoe  is  assumed  to  be  rigid. 

d.  The  friction  coefficient  is  a  linear  function  of 
temperature  and  it  does  not  vary  with  pressure,  wear 
and  anvironment. 

2.  Pressure  Concept 

Td  analyze  an  internal  shoe  refer  to  Fig.  2,  which 
shows  a  shoe  pivotal  at  a  fixed  point  with  the  actuating 
force  acting  at  the  other  end  of  the  shoe.  The  mechanical 
arrangement  does  not  permit  pressure  to  be  applied  at  the 
pivot,  therefore  the  pressure  at  this  point  is  zero.  If  the 
shoe  rotates  through  a  small  angle  about  A,  the  radial 
movement  of  any  point  on  the  arc  of  contact,  is  proportional 
to  the  moment  arm  of  this  point  from  the  pivot.  Assuming 
that  the  material  of  the  crake  lining  and  support  obey 
Hooke's  law,  the  pressure  at  this  point  will  also  be 
proportional   to   this   moment   arm.    The    distance   is 
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proportional  to  sine  .  Therefore,  the  relations 
pressure  at  any  point  and  the  maximum  pressure,  p  , 
given    by    the    following    formula; 

P  Po 


sini 


between 
will  be 


(1) 


a 

From  this  formula  it  can  be  seen  that  the  frictional 
material  at  the  heel,  contributes  very  little  to  the  braking 
action,  therefore  it  is  better  to  begin  the  friction 
material  at  an  angla  0!  greater  than,  say  0.15  rai.  It  can 
be  seen  also  that  the  pressure  will  be  maximum  when  6  =90° 
or  if  the  toe  angle  92  is  less  than  90°,  then  the  pressure 
will  be  naximura  at  the  toe.  For  good  performance  it  is 
recommended  to  concentrate  as  much  frictional  material  as 
possible  in  the  neighborhood  of  the  point  of  maximum 
pressure   [Ref.  2]. 

3.   Actuating  Force  and  Torgue  Calculation 

From  Fig.   2,   it  can   be  seen  that  the  differential 
normal  force  on  an  elament  area  of  the  lining  will  be; 


dN  =  pdA 


(2) 


where  dA  is  an  area  element  of  tha  lining  and  it's  magnitude 
is; 


dA  =  rbdi 


(3) 


In  Equation  3,  r  is  the  inside  drum  radius  and  b  is  the  drum 
width.   Substituting  for  p  and  dA  gives; 


p  brsin9 
dN  =   a  .  s —  d9 


(4) 


At  the  same  point  the  differential  frictional  force  is; 

df  =  ydN  (5) 

where  u   is  the  coefficient  of  friction. 
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rhe  actuating  force,  F,  can  be  calculated  using  the  fact 
that  the  summation  of  the  moments  about  the  hinge  pin  is 
zero.      The    moment   due    to   frictional    forces    is; 

M„    =fl;02(r-cos6)df  (6) 

i      y  i 

where  a  is  the  distance  from  the  pivot  to  the  center  of 
rotation.  Substituting    the      value    of      dF    and      integrating 

from   9j    to    e2    gives; 

yp    br 2 
Mf    =  jj^ — {(cosei-cose2)+|?(sin26i-sin262)}  (7) 

a 

where  y      is    assumed    to    be    constant    along   the    lining. 

Similarly    the   moment    due    to    normal    forces   is    given    by; 

M      =    /92asin8dN  (8) 

n      y  i 

Substituting  the        value  of  dN        and  integrating 

from    G i    to    9  2    gives; 

p  bra 
Mn  =  sin9  {0-5(9 2"9 i)-0-25( sin9 2-sin6i)  (9) 

a 

The  actuating  force  must  balance  the  moments, therefore; 

p  =  _2 I 

d  (10) 

where  d  is  the  distance  from  the  hinge  to  the  point  of 
application  of  F.  The  torque  applied  to  the  drum  by  the 
brake  shoe  is; 

To  =e/92rdf  (11) 

After  substituting  the  value  of  df  and  integrating  ; 

yp  br2 
T0  =  — t-^-k — (cos6i-cos92)  (12) 
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H.  Rate  of  Heat  Generated  ani  Deceleration  Calculation 
The  differential  rate  of  heat  generated  by  an 
element  area  of  the  lining  is  equal  to  the  velocity  of  the 
inside  surface  of  the  drum  relative  to  the  lining  ,  rimes 
the  differential  frictional  force  acting  on  the  element 
area; 


dQ  =  Vrdf 


(13) 


Assuming   the    brake   is    on   a    vehicle    wheel   with    a   radius   of   R, 
the    inside      surface   velocity   is   equal    to; 


V        =    ±rV 

r        R 


(14) 


where  V  is  the  velocity  of  the   vehicle  and  is  a  function  of 

time. 

If    7  =  7  (t)    then    7    =7    (t)       and   the    heat    generated    will   be    also 
r      r 

a   functioa    of  time.    Substituting    the    values   of   7      and   df   and 

3  r 

integrating  from  Q1    to  e2  *    *e  9et  the  following  formula  for 
the  heat  generated  at  any  time  t. 


p  by    2 
Q(t)  =  Sin9  (|0(cos9i~cos92)V(t) 
a 


(15) 


The  kinetic  energy  of  a  vehicle  of  weight  H  is  given  by; 

E  =  M)V2  (16) 

Z   6 

Mote  that  if  the  brake  is  on  a  four  wheel  vehicle,  there 
will  be  eight  shoes.  Assuming  all  are  leading  shoes,  each 
will  stop  one-eight  of  the  vehicle  weight,  so  W/8  must  be 
used  in  Equation  (16)  .  The  rate  of  change  in  the  kinetic 
energy   is; 


dE 


W,„dV 


S  -  (i)vdt 


(17) 


From  the  energy   conservation  law  the  rate  of   change  in  the 
kinetic  energy  is  equal  to  the  heat  generated; 
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!(t)  = 


dE 
dt 


(18) 


Substituting  the  value  of  Q  (t)  and  dE/dt,  it  is  seen  that 
the  velocity  7(t)  cancels  and  so  the  deceleration  is  not  a 
function  of  time.   Therefore  the  deceleration,  dc,  is; 


dc  =  dt 


Cf>ifS-<|)(c^e»-cose*) 


(19) 


The  velocity  at  any  time  is; 


V  =  V.-dct 

1 


(20) 


where  7.  is  the  initial  velocity.  Substituting  the  velocity 
ii  Equation  (16) ,  yields  the  rate  of  heat  generated  as  a 
function  Df  time, 


paby   r2 
Q(t)  =  _|_-(|)(cos0i-cos82)(Vi-dct) 


(21) 


In  this  study  the  friction  coefficient  was  taken  as  constant 

up   to  a   temperature  of   90°C   and  after   90°C,    decreases 

linearly  to  zero  at  a  specified  temperature,  Tmr: 

max  . 

,  u  T  <  90  °C 


AT 


(T-90) 


I  0 


90  °C  <  T  <  T 

—  max . 


T  >  T 


max. 

where  y   is  the  cold  coefficient  of  friction  and  y] 
hot  coefficient  of  friction. 


(22) 


is  the 


E.   SURFACE  TEMPERATURE  CALCULATION 

Since  the  function  of  a  brake  is  to  convert  kinetic 
energy  into  heat,  surface  temperatures  of  brake  linings  and 
drums  are  lost  important.  Therefore  it  is  necessary  to  know 
the  temperature  of  the  mechanism  during  and  after  any  stop. 
The  temperatures  were  calculated  by  the  finite  difference 
method. 


1.  Assumptions 

a.  One  dimensional  heat  flow-Ihe  heat  flow  is  from  the 
inner  surface  to  the  outer  surface  of  the  drum. 

b.  Constant  heat  transfer  coefficient. 

c.  No  heat  dissipated  by  radiation. 

d.  The  heat  is  generated  on  the  inner  surface. 

2.  Temperature  Analysis 
a.   Theory 

The  differential  equation  to  be  solved  in  order 
to  find  the  temperature  in  the  drum,  based  on  the 
assumptions,  is; 

ifr  +  Q  =  (i)|T  (23) 

with   the   following   boundary   conditions: 

at   x=0         heat    is    generated, 

at   x=t]c        heat      is      transferad      to      the      atmosphere      by 

convection. 

In  the  equation   above  Jc  is  the   thermal  conductivity,  a  is 

the   thermal  diffusivity,    t  is   time   and  tk   is  the   drum 

thickness.     This  ejuation   can   be   solved  by   the   finite 

difference  method   "Ref.  3].    The   finite  difference   model 

used  here  is  shown  in  Fig.  3.    The  rate  of  change  with  time 

of  the  internal  energy  of  a  node  i  is  approximated  by; 
TP+1_TP 

||  =  PcAV„-i^  (24) 

where  p  is  the  density,   c  is  the  specific  heat  and  VQ   is 

the  drum  volume. 

Now  define  the  thermal  capacity  as 

C.  =  p.c.AVo-  (25) 

l    "i  l    l 

The   forward    difference      equation    for   all   nodes      and   boundary 
conditions    is: 


19 


rpP_rpi 


P+1_TP 


where    B..  is   the    thermal    resistance 

th,ij  +1 

Solving   the    above   equation    for   T"  gives 


(26) 


JP+1   - 


th,ij      1  i      th.ij 


(27) 


The  thermal    resistanca      can    be  calculated   from      the    geometry 
and  boundary  conditions  [Ref.    3].  To  ensure 

stability      At      must      be    equal      or      less      than    the      following 

nodal  relation; 

C. 


At    <    ( 


I- 


(28) 


th,ij 
With  the   assumptions  made,    the  drum  can   be  viewed   as  an 

infinite  plate,   with   heat  generated  at  the   surface  of  the 

first  node,  as  shown  in  Pig.  3.   It  is  assumed  that  in  every 

drum,   there  are  two  shoes  and   that  both  are  leading  shoes. 

Therefore,  two  times  Q.  must  be  taken. 


,P+1 


TP 


i    =  (2QP+£ 


It 


R.  .   .  .  C. 
th,ij   i 


+  (1-  gfcri-) 


(29) 


i  " th, ij 
b.   Formulation 

In  the  computer  program  5  nodes  were  taken.  In 
order  to  check  accuracy,  the  program  was  run  with  7  and  10 
nodes.  In  each  case  the  result  was  the  same  within  5  C. 
The  heat  is  generated  in  the  inner  drum  surface.  Therefore  Q 
appears  in  the  formula  of  temperature  in  the  first  node  and 
for  all  the  other  nodes  Q  is  equal  zero.  With  the 
assumptions  mentioned  above,  the  heat  transfer  through  the 
drum  is  solved  as  a  heat  transfer  problem  through  an 
infinite  plate,  with  heat  generation  at  the  inner  surface 
and  with  a  heat  convection  boundary  on  the  outer  surface  as 
shown  in  Fig.  3.  Equation  (29)  can  be  simplified  using  two 
dimensionless  parameters,  Biot  and  Fourier  modulii, 
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B.  =  ^  (30) 

1     k 

F0  =  aAt  (31) 

(Ax)z 

The  final  aquations  for  calculating   the  temperatures  at  the 
nodes  now  become; 
For  tha  first  noda; 


2QPAt 


P+1  =  -£±—  +  (1-2F0)TP  +  2F0TP 


1 

For  tha  interior  nodas; 


(32) 


TP+1  =  F0{TP  -  +  T?+1  +  (i-  "  2)TP}                (33) 

1             1-1      1  +  1  r  o         1 

For  tha  last  node; 

TP+1  =  2F0(TP  -+B.T  +  (~-  -B.-DTP}               (34) 

n           n— l   i  oo  zr  o    i     n 


F.   BRAKE  DUTY  CYCL2 

In  addition  to  tha  parameters  mentioned  above  the  design 
of  a  braka  depends  on  the  initial  speed,  final  speed,  number 
of  stops,  and  the  rest  time  between  each  stop.  In  this 
analysis  a  general  duty  cycle  was  considered  so  that  the 
initial  speed,  final  speed  ani  the  acceleration  period 
between  stops  can  be  different  for  each  part  of  the  design. 

In  tha  design  examples  presented  here,  a  vehicle  was 
stopped  four  consecutive  times  with  the  following  cycle; 

Rest 


sec. 

20.0 
20.0 
20.0 


Initial 

Final 

Speed 

Speed 

m/sec. 

m/sec 

1 

25.0 

0.0 

2 

25.0 

0.0 

3 

25.0 

0.0 

4 

25.0 

0.0 
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III.   OPTIMIZATION 

A.   INTRODUCTION 

Engineering  analysis  using  the  digital  computer  has 
become  commonplace.  It  is  less  common  to  use  the  computer 
to  make  the  actual  design  decisions,  such  as  sizing  of 
structural  members  or  placement  of  mechanical  linkages.  This 
may  be  largely  attributed  to  the  fact  that  fully  automated 
design  requires  techniques  that  are  unfamiliar  to  much  of 
the  engineering  community. 

In  many  enginaering  problems,  it  is  necessary  to 
determine  the  minimum  or  maximum  of  a  function  of  several 
variables,  limited  by  various  linear  and  nonlinear 
inequality  constraints.  It  is  seldom  possible,  in  practical 
applications,  to  solve  these  problems  directly,  and 
iterative  methods  are  used  to  obtain  the  numerical  solution. 
Machine  calculation  of  this  solution  is,  of  course, 
dasirable.  The  CONMIN  program  is  available  to  solve  a  wide 
variety  of  such  problems  [Ref.  4]. 

CONMIN  is  a  FORTRAN  program,  in  subroutine  form,  for  the 
minimization  of  a  multi-variable  function  subject  to  a  set 
of  inequality  constraints.  The  basic  optimization  algorithm 
is  the  Method  of  Feasible  Directions  [Ref.  5].  The  user  must 
provide  a  main  calling  program  and  an  external  routine  to 
evaluate  the  objective  and  constraint  functions  and  to 
provide  gradient  information.  If  analytic  gradients  of  the 
objective  or  constraint  functions  are  not  available,  this 
information  is  calculated  by  finite  difference.  While  the 
program  is  intended  primarily  for  efficient  solution  of 
constrained  problems,  unconstrained  function  minimization 
problems  nay  also  be  solved,  and  the  Conjugate  Direction 
Method  of  Fletcher  and  Reeves  is  used  for  this  purpose 
[Ref.  6]. 
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B.   DEFINITION  OF  TERMS 

Most  lisciplines  have  a  unique  set  of  nomenclature  used 
to  describe  the  concepts  within  that  discipline.  Some  of  the 
commonly  used  terms  in  numerical  optimization  are  summarized 
hare. 

Objective-  The  value  of  the  function  which  is  to  be 
minimized  or  maximized  during  the  optimization  process. 
Synonyms  are  cost,  merit  and  payoff.  The  common  mathematical 
designation  is  F  (X) .  In  the  present  study  the  objective  was 
to  minimize  the  material  in  the  brake  drum. 

Design  variables-  Th*  parameters  to  be  changed 
during  the  optimization  process  in  order  to  minimize  or 
maximize  the  value  of  the  objective  function.  Synonym; 
decision  variables.  The  common  mathematical  designation  is 
the  vector  X.  Design  variables  considered  in  this  study 
include,  drum  thickness,  width,  the  angle  between  the  hinged 
pin  and  the  end  of  the  lining,  and  the  distance  from  the 
pivot  to  the  center  of  rotation. 

Inequality  constraints-  Dne-sided  conditions  which 
must  be  mathematically  satisfied  for  the  design  to  be 
acceptable.  The  common  mathematical  term  is  G(X)<0  or 
G(X)>0.  If  the  ineqality  condition  is  satisfied  on  G ( X)  ,  the 
design  is  acceptable,  (feasible) .  If  it  is  not  satisfied, 
the  design  is  not  acceptable  (infeasible) .  Constraints 
considered  here  include,  vehicle  stopping  time,  maximum  drum 
temperature,  and  actuating  force. 

Side  constraints-  Upper  and  lower  bounds  on  the 
individual   design  variables   X  .    The  common   mathematical 

representation  is  X.<  X.<  X.. 

c  111 

Design    space-    The      n-dimensional   mathematical   space 
spanned   by    the   vector    of    design    variables   X. 

Active  constraint-  Constraint  G . (X)  is  called  active 
if  its  value  is  zero  (or  near  zero  for  computational 
purposes)  . 
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G  (X)<0, 

J 

3  (X)  >0. 

J 


Inactive  constraint-  Constraint  G. (X)  is  inactive  if 

J 

Violated  constraint-  Constraint  G.(X)  is  violated  if 

J 


C.   THE  OPTIMIZATION  PROCESS 

The  general  design  optimization  problem  can  be  stated 
mathematically  as  follows:  Find  the  set  of  variables  X., 
i=1,2....a,  which  will 

Minimize  F  (X  ) 
Subject  to: 


(35) 


G.(X)  <   0     j-1,2 m 

X1  £  X.  <  XU     i=l,2 n 


(36) 

(37) 


Vector  X  contains  the  set  of  independent  design 
variables  X.,  1=1,2. ...n.  X  may  represent,  for  example 
width,  thickness,  and  angles  in  the  brake  optimization.  The 
objective  function  used  here  is  the  drum  volume. 

Equation  (36)  defines  the  in=qality  constraints  imposed 
on  the  design.  For  example,  if  the  temperature  on  the  inner 
drum  surface  must  not  exceed  a  specified  value  I,  the 
associetei  design  constraint  becomes,  in  normalized  form 


-—■  -  1  <  0 


(38) 


The  lDwer  and  upper  bounds  on  the  design  variables, 
given  by  Eg.  (37)  ,  limit  the  region  over  which  the  functions 
F(X),  and  G  (X)  are  defined.  These  constraints  are  often 
referred  to  as  side  constraints  because  they  form  the  sides 
or  bounds  of  the  n-dimensional  spaced  spanned  by  the  design 
variables  X. 

If  all  the  inequalities  of  Sqns.  (36)  and  (37)  are 
satisfied,  the  design  is  said  to  be  feasible;  if  any  of 
these   CDnditions  are   not   satisfied,    the  design   is   not 
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feasible.  If  F(X)  is  a  minimum  and  the  design  is  feasible, 
it  is  also  optimum,  or  at  least,  a  relative  optimum.  Note 
that  because  the  objective  and  constraints  may  be  nonlinear, 
there  may  be  multiple  minima  in  the  design  space  that  cannot 
be  identified  using  current  methods.  While  this  is  a  matter 
for  concern,  since  it  is  desired  to  find  the  true  optimum, 
it  must  be  remembered  that  the  same  mathematical  conditions 
exist  if  the  design  process  is  not  automated.  However,  using 
optimization  techniques,  it  is  a  simple  matter  to  restart 
the  optimization  from  several  initial  points  in  the  design 
space  and  thereby  improve  the  probability  of  obtaining  the 
tcue  optiaum  design,  a  process  that  would  be  quite  time- 
consuming  in  manual  design. 

Equations  (35) -(37)  define  the  nonlinear  constrained 
optimization  problem.  If  Egs. (36}  and  (37)  are  not  imposed 
on  the  design,  the  optimization  problem  is  defined  by 
Eq.(35)  alone  and  is  therefore  an  unconstrained  minimization 
problem. 

Host  nonlinear  optimization  algorithms  update  the  vector 
of  design  variables  by  the  iterative  relationship; 


r 


q-l  + 


aS 


5q 


(39) 


where  g  is  the  iteration  number,  vector  S  is  the  direction 
of  search  in  the  design  space,  ani  the  scalar  a  is  referred 
to  as  a  aove  parameter  which, together  with  S,  determines  how 
much  the  vector  X  is  changed  during  the  q-th  iteration.  An 
initial  design  defined  by  X  must  be  supplied.  The 
optimization  process  then  proceeds  in  two  steps.  First, 
the  direction  S,  which  improves  the  design,  is  found,  and 
second,  the  scalar  a  ,  is  determined  which  improves  the 
design  as  much  as  possible  when  moving  in  this  direction. 
The  process  is  repeated  until  there  is  no  further  design 
improvement,   indicating  that  this  is  the  optimum  attainable 
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design.  Por  farther  details  see  Ref.  7. 

D.   COPES  AND  SUBROUTINE  ANALIZ 

In  order  to  simplify  the  use  of  CONMIN  and  to  further 
aid  in  the  design  optimization  process  a  Control  Program  For 
Engineering  Synthesis,  COPES,  was  developed  by  Vanderplaats 
[Ref.  7],  COPES  is  the  main  program  (recall  that  CONMIN  is 
written  in  subroutine  form) .  The  user  must  supply  an 
analysis  subroutine  with  the  name  ANALIZ,  which  will 
calculate  the  various  parameters.  This  subroutine  has  three 
segments;  INPUT,  EXECUTION,  OUTPUT. 

All  parameters  which  may  be  design  variables,  objective 
functions  or  constraints  are  contained  in  a  single  labeled 
common  block  called  GLOBCH. 

Copes  Terminology 

The  ZOPSS  program  currently  provides  six  specific 
capabilities; 

1.  Simple  analysis,  just  as  if  ZOPES  was  not  used. 

2.  Optimization-Minimization  or  maximization  of  one 
calculated  function  with  limits  imposed  on  other 
functions. 

3.  Sensitivity  analysis-  The  effect  of  changing  one  or 
more  design  variables  on  one  or  more  calculated 
functions. 

4.  Two-yariable  function  space-Analysis  for  all  specified 
combinations  of  two  design  variables. 

5.  Optinum  sensitivity-  The  same  as  sensitivity  analysis 
except  that,  at  each  step,  the  design  is  optimized  with 
respect  to  the  independent  design  variables. 

6.  Approximate  optimization-  Optimization  using 
approximation  techniques.  Usaally  more  efficient  than 
standard  optimization  for  up  to  10  design  variables  or 
if  multiple  optimizations  are  to  be  performed  [Ref.  7]. 
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17.       DESCRIPTION    OF    TffE    COMPUTER    PROGRAM 

A.  GENERAL  PROGRAM  ORGANIZATION 

A  functional  block  diagram  of  the  program  is  presented 
in  Fig.  H .  A  general  description  of  the  subroutines 
contained  in  the  program  is  given  here.  Appendices  A 
through  D  discuss  the  preparation  of  input  data,  list  the 
important  computer  program  nomenclature,  and  list  the 
program. 

B.  SUBROUTINES 

1.  Subroutine  ANALIZ 

Subroutine  ANALIZ  organizes  the  basic  analysis  used 
in  the  optimization.  It  controls  the  reading  of  the  initial 
design  description  and  calculation  of  the  values  of  the 
objective  function,  constraints,  and  all  other  parameters 
necessary  to  solve  the  problem.  COPES/CONMIN  updates  the 
design  to  minimize/maximize  the  objective  function, 
iterating  until  no  further  improvement  in  the  objective 
function  is  possible  without  violating  one  of  the 
constraints.  COPES/CONMIN  calls  subroutine  ANALIZ  to  obtain 
the  function  value  iucing  the  optimization. 

2.  Subroutine  INPUT 

This  subroutine  reads  all  input  data  associated  with 
the  brake  analysis.  Instructions  for  problem  deck 
preparation  are  given  in  appendix  B. 

3.  Subroutine  TEMPR 

This  subroutine  calculates  the  heat  transfer 
constants  such  as  the  thermal  capacity  of  each  node  and  the 
resistance  of  each  node,  determines  the  time  increment  in 
order  to  insure  a  stable  solution,  and  calculates  the  rate 
of  heat  generation.  In  order  to  calculate  the  temperature  of 
each  node,  it  calls  two  subroutines.  From  subroutine  BRAK  it 
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obtains  the  deceleration  needed  to  calculate  the  rate  of 
heat  generated  and  from  subroutine  TEMA  it  obtains  the 
temperature  rise  of  each  node.  Then  it  calculates  the 
temperatures  during  the  time  that  the  brake  is  not  in  use. 
This  subroutine  is  also  capable  of  calculating  the 
temperature  rise  of  a  drum  when  a  constant  rate  of  heat 
dissipation  is  given. 

4.  Subroutine    TE3A 

This  subroutine  calculates  the  temperature  of  each  node. 
As  mentioned  before,  tne  heat  is  generated  on  the  inner 
surface,  and  on  the  outer  side  of  the  drum  the  heat  is 
dissipated  by  convection.  The  formulas  used  were  developed 
by  the  finite  difference  method,  and  are  given  in  section 
II-E-2. 

5.  Subroutine  BRAK 

This  subroutine  calculates  the  torque,  actuating  force, 
and  the  friction  moment  of  one  shoe.  It  also  calculates  the 
drum  volume  and  the  deceleration  of  the  machine.  The 
subroutine  takes  into  consideration  a  constant  friction 
coefficient  until  a  temperature  of  90°C  is  reached  and  a 
linear  decrease  in  the  friction  coefficient  for  higher 
temperatures.  More  details  are  given  in  section  II-D-4. 

6.  Subroutine  DUTPUT 

This  subroutine  echos  the  input  data  and  prints  out 
tne  thermal  and  mechanical  information  for  the  brake.  An 
example  of  the  output  obtained  from  this  subroutine  is  shown 
in  Table  1  and  Table  2. 


V.  T5ST  PROBLEM  AMD  RESULTS 

The  computer  program  was  tested  with  the  data  specified 
in  Table  t.  The  objective  function  which  was  minimized  was 
the  volums  of  the  drum  material.  Design  variables  were  the 
drum  width,  the  angls  between  tha  hinged  pin  and  the  end  of 
the  lining,  the  ratio  of  the  pivot  to  center  of  rotation 
distance  to  drum  radius,  and  the  drum  thickness.  The  side 
constraints  (limits)  on  the  design  variables  were; 

Upper 
Bound 


Design 

Lower 

Variable 

Bound 

1. 

(3}     Width, b 

0.0 

2. 

(5)     Theta    2, 

1.2    rad. 

3. 

(12) Ratio, Rd 

0.1 

4. 

(1B)  Thickness,  tk 

40      mm. 

80  mm. 

2.5  rad. 

0.  9 

No  bound 

The  number   in   parentheses   is   the   location   of   the 

variable  in  the  COMMON  block  in  the  computer  program. 

Constraints  were  imposed  on  the   actuating  force  F,   the 

maximum  temperatura,  T     ,  on  the  inner  surface  of  the  drum 
r  max 

and  stopping  time,   t. 


Constrained 

Lower 

Variable 

Bound 

1. 

(9)     Force, F 

200.0    N- 

2. 

(25)     Time,    t 

No    bound 

3. 

('4)     Temperaturs 

*T 

No    bound 

Upper 

Bound 

2500.0    N-m 

7.00 

Sec. 

230.0 

°C 

The  vehicle  which  weights  25700.0  Newtons  is  stopped 
four  consecutive  times  from  a  velocity  of  90.0  Km/hr  to 
zero,  with  an  acceleration  period  of  20.0  sec.  between 
stops.  The  values  of  the  design  variables  and  the 
constraints  before  and  after  optimization  are; 
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Objective 

Function 

Drum  Volume 

Design 

variables 

Width 

Theta  2 

a/r 

Thickness 

Constraints 


Before 

3pt imization 


0.754    E-03    m 


0.08    a 

2.  10    Rad, 
0.75 
0.010    a 


Actuating 

Force  2815.1    N 

Stopping 

time    (last    stop)    7.04    sec. 

Temperature 

after   last    stop      348. 7 °C 


After 
Optimization 


0.159    S-02    m 


0.08    m 
1.92    Rad. 
0.755 
0.020    m 


2086.3    N 


7.00    sec. 


229. 2°C 


Note  that  the  objective  function  increased  as  a  result 
of  optimization.  This  is  because  the  initial  design  violated 
constraints    on   stopping   time   and   maximum   temperature. 

Further  results  are  listed  in  Tables  1  and  2.  In 
addition  to  optimization,  a  sensitivity  analysis  of  the 
design  variables  and  a  two-variable  function  space  analysis 
for  width  and  thickness  were  performed.  The  graphical 
results  are  given  in  Figs.  5  through  17.  The  results  can  be 
summerized    as   follows; 

a.    The   effect    of   changing   the      inside  drum   radius   with   all 
other    design   variables   held   constant; 

As  shown  in  Figs.  5-7,  for  small  inside  drum  radii  the 
drum  temperature  is  very  high.  The  stopping  time  is 
long  and  the  torgue  is  low.  Inside  drum  radii  over  130 
mm  give  reasonaole  drum  temperature  and  stopping  time, 
for    the    example    considered. 
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b.  As  seen  in  Figs.  8-10,  the  affect  of  changing  the  drum 
width  with  all  other  design  variables  held  constant  is 
the  same  as  described  above. 

c.  The  affect  of  changing  tha  drum  thickness  with  all 
other  parameters  held  constant  is;  For  a  drum  thickness 
up  to  6  mm,  the  stopping  tiae  and  drum  temperature  are 
consiierably  high.  Over  16  urn  thickness,  the  stopping 
tima  re  nains  almost  constant.  For  a  small  thickness  the 
torqua  is  very  low  due  to  the  high  temperatures.  For 
thicknesses  over  20  mm,  the  torque  remains  about 
constant. 

d.  The  effect  of  changing  the  angle  between  the  hinged  pin 
and  the  end  of  the  lining  is;  For  a  small  92  angle  the 
stopping  time  is  very  long  because  the  torque  is  low. 
The  stopping  tima  becomes  reasonable  when  02  >1.8  Rad. 
DbviDusly  there  is  an  increase  in  the  drum  temperature 
as  G2  increases  but  the  overall  change  in  temperature 
is  small. 

e.  From  the  two  variable  function  space.  Fig.  17,  it  can 
be  seen  that  the  constant  volume  line  and  the  constant 
temparature  line  are  almost  parallel,  this  leads  to  the 
conclusion  that  for  the  cycle  taken,  the  drum  is  a  heat 
sink,  and  the  amount  of  heat  dissipated  by  convection 
during  this  cycle  is  small. 
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VI.   TEHPE5ATUBE  RAISE  -  SIMPLIFIED  CALCULATION 

A  simplified  way  of  finding  the  temperature  rise  of  the 
drum  is  by  using  the  equation; 

Q  =  |cAT  (40) 

and  setting  Q  equal  to  the  amount  of  heat  generated  using 
Equation  (21)  from  section  II-D-4.  This  equation  is  in 
common  use  in  engineering  texts  (See,  for  example,  Refs.  2 
and  8) .  The  temperature  rise  calculated  this  way,  is  the 
average  temperature  of  the  drum,  and  not  the  temperature  on 
the  interface,  which  can  be  much  higher  (depending  on  the 
rate  of  heat  generated) .  Extreme  temperature  gradients  cause 
distortion  and  excessive  surface  wear.  Therefore  it  isn't 
always  acceptable  to  use  the  simplified  formula.  From 
experience,  it  has  been  found  that  the  surface  wear 
increases  dramatically  as  interface  temperatures  approach 
400  to  500°F  (205  to  260  °C),  [Ref.  8]. 

A  comparison  of  the  temperatures  calculated  on  the  inner 
drum  surface,  outer  drum  surface,  and  the  average  drum 
temperature  calculated,  using  equation  (40) ,  is  given  in 
Fig.  18.  The  graph  shows  the  temperature  rise  for  a  vehicle 
stopped  from  a  velocity  of  90  sm/hr.  From  this  graph,  it  can 
be  seen  that  the  drum  temperaturs,  based  on  equation  (40), 
after  the  vehicle  stopped  is  about  the  average  temperature 
of  the  inner  and  outer  surface  temperatures. 

The  results  show  that  the  drum  will  reach  an  uniform 
temperature  of  about  53°C,  in  15  sec,  after  the  vehicle  has 
stopped. 

Calculating  the  temperature  with  the  simplified  formula, 
can  lead  to  errors  in  the  time  needed  to  stop  the  vehicle. 
Because  the   temperature   calculated   with   the   simplified 
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formula  is  lower  than  the  temperature  at  the  friction 
interface,  the  calculated  friction  coefficient  is  higher 
than  the  actual  friction  coefficient.  Therefore  the 
calculate!  stopping  time  will  be  shorter  than  the  real 
stopping  time.  Ail  this  is  true,  provided  the  friction 
material  behaves  as  assumed  in  section  I-D-4. 

Because  high  temperature  is  detrimental  to  both  the 
stopping  ability  and  the  wear  characteristics  of  the  brake, 
it  is  important  that  the  interface  temperature  be  calculated 
with  reasonable  accuracy  in  design.  Fig.  18  clearly  shows 
the  temperature  differences  resulting  from  the  two 
approaches. 

This  difference  in  results  is  compounded  when  the 
simplified  equation  is  used  for  design.  Table  3  presents  the 
design  results  based  on  the  simplified  approach.  This  design 
represents  an  apparent  material  savings  of  27%.  However, 
when  this  optimum  is  analysed  using  the  finite  difference 
heat  transfer  solution,  the  maximum  temperature  is  268.8 °C 
and  the  Last  stopping  time  is  7.54  sec.  This  time  violates 
the  constraint  by  7.7%,  Perhaps  more  importantly,  the 
temperature  at  the  interface  of  about  269  °C  would  surely 
lead  to  premature  failure.  Therefore  this  design  is  clearly 
too  unconser vative  to  be  acceptable. 
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711.  CONCLUSION 


In   summary,   a   numerical   optimization  program  is   an 

effective   way   of   finding  a   solution   to  an   engineering 

problem,  provided  reasonable  care  is  used  in  formulating  the 
problem. 


VIII.  RECOMMENDATIONS  FOR  FUTURE  INVESTIGATIONS 

The  study  has  shown  the  feasibility  of  using  numerical 
optimization  in  the  design  of  Internal-Expanding  Rim  Brakes 
with  two  Leading  shoes.  Further  studies  on  the  same  design 
may  be  pursued  by  eliminating  some  of  the  restrictions.  For 
example; 

1.  To  add  heat  dissipation  by  radiation. 

2.  To  investigate  drum  temperatures  for  a  drum  with  fins. 

3.  To  take  into  consideration  changes  in  the  surface 
pressure  as  a  function  of  friction  coefficient. 

4.  To  repeat  all  the  calculations  for  a  drum  in  which 
there  is  one  trailing  shoe  and  one  leading  shoe. 

5.  To  add  the  effect  of  centrifugal  forces  for  clutches. 
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APPENDIX  A 


LIST  OF  PARAMETERS 


A  complete  listing  and  description  of  all  variables  used 
ia  the  prsgram,  is  not  practical.  The  variables  listed  in 
this  appendix  are  common  to  several  subroutines  of  the 
program  and  will  assist  the  reader  in  a  study  of  the 
program.  The  Global  location  is  the  location  of  the 
parameter  in  the  common  block  called  GLOBCM.  This  common 
block  is  the  means  by  which  information  is  transfered 
between  the  subroutines  and  the  CDPES/CONMIN  program. 


Global      Fortran    Math.      Definition 
Location     Name      Symbol    


1  HI        r         Inside  drum  radius   (m) 

2  T         TQ         Torque  of  one  shoe  (N-m) 

3  WDTH       b  Drum  width   (m) 

U  PRSA      p         Pressure  between  lining 

and  drum   (N/ra2) 

5  TETA1      0i         The  angle  between  the 

hinged  pin  and  the  (Rad.) 
begining  of  the  lining 

6  TETA2     0  2         The  angle  between  the 

hinged  pin  and  the  end 
of  the  lining   (Rad.) 

7  FRMNT     Mf        Friction  moment   (N-m) 

8  AN13RT     Mn        Normal  momemt   (N-m) 

9  ACFRC     F         Actuating  force   (N) 

10  C         d         Distance  from  actuating 

force  to  the  hinged  pin(m) 
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11  Q         3  Heat  generated   (J/sec.) 

12  RD        a         Distance  from  pivot  to 

center  of  rotation   (m) 

13  CMIU      u         Cold  friction  coefficient 

c 

14  HMIU  u  Hot    friction   coefficient 

h 

15  AHIU      u  Friction  coefficient  at 

any  temperature 

16  SRFC  Drum  surface  area  (m2) 

17  RO  Outside  drum  radius  (a) 

18  THK       tk        Drum  thickness  (m) 

19  DX  An  incremental  thickness 

(m) 

20  RUES  R  Wheel   radius    (m) 

21  W  W  Car's   weight    (N) 

22  DCCE  dc  Deceleration       (m/sec.?  ) 


Q 

Q 

RD 

a 

cmiu 

HMIU 

u 

AMIU 

SRFC 

RO 

THK 

tk 

DX 

RHER 

R 

W 

VI 

DCCE 

dc 

TOT 

ECEN 

a 

NWRT 

TIME 

t 

VOL 

Vo 

TEPL(5) 

T 

HWR 

NWRA 

NWRQ 

NEL 

NSEG 

?I 

IT 

G 

g 

K 

k 

23  TOT  Total   time       (sec.) 

24  ECEN  a  Ecentricity    (m) 

25  NWRT  Write   statement  control 

26  TIME  t  Time    (sec.) 

27  70L  Vo  Drum    volume    (m3) 
28-32                 TEPL(5)         T                      Temperature   at   time   p+1 

(sec.) 

33  NWR  Write   statement   control 

34  NWRA  Write  statement   control 

35  NWRQ  Write   statement   control 

36  NEL  Number   of    elements 

37  NSEG  Number   of    segments 

38  ?I  it  Constant 

39  G  g  Gravitational   constant 

40  K  k  Thermal   conductivity 

(J/m-C) 

41  HCNV      h         Convection  heat  coef f . 

(W/m2-C)  . 

42  SPHT      c         Specific  heat  (J/Kg-°C) 
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43 

RHO 

P 

44 

DTAO 

45 

DFTM 

T 

46-50 

TEMP 

T 

51-57 

RES 

58-63 

TC 

C 

64 

BIO 

3i 

65 

FUR 

Fo 

66-72 

MVT 

73-79 

VT 

80 

tfELO 

81 

NELT 

82 

TETAA 

e 

83 

84 
85 
86 


ACOF 
TISI 
ZMAN 
NS  HO 


Density  (Kg./m3) 
Time  increment  (sec.) 
flax.  temp,  difference ( °C) 
Temp,  at  time  p  (sec.) 
Heat  resistance  ( °C/J) 
Heat  capacity   (J/°C) 
Biot  moduli 
Fourier  moduli 
Control  parameter 
Control  parameter 
Number  of  elements+1 
Number  of  element s+ 2 
The  angle  at  which  the 
pressure  between  the 
lining  and  drum  is 
maximum.  (Rad.) 
Constant 

Initial  temperature ( °C) 
Time  increment  (sec.) 
Number  of  shoes 
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APPENDIX  B 
INSTRUCTIONS  FOR  PROBLEM  DATA  PREPARATION 

Although  the  procedure  is  straight  forward,  preparation 
of  input  data  for  the  program  requires  attention.  Errors  are 
easy  to  make  and  difficult  to  locate.  Input  data  is 
described  here  for  the  brake  analysis.  For  instructions  on 
data  preparation  for  optimization  see  Ref.  7.  Input  data 
should,  in  general,  follow  the  steps  outlined  below.  The  use 
of  the  standard  FORTRAN  Eighty  Columm  Coding  Sheet  is 
recommended.  Integer  constants  aust  be  right  justified  in 
the  appropriate  field.  There  are  eight  input  cards,  read  by 
subroutine  INPUT,  to  describe  the  initial  design,  material 
properties  and  constants.  Card  format  is  given  in 
parenthesis  followed  by  specific  instructions  where 
necessary. 

1.  First  Card  (110)  -  Duty  cycle  information. 

Cols  1-10   :  Total  number  of  consecutives  stops 
and  accelerations  (NSEG) 

2.  Second  Card  (I10,3F10.0)  -  Duty  cycle  information, 
a.   Cols  1-10   :  Control  number. 

1  means-deceleration, 

2  means-brake  not  in  use, 
Velocity  at  start  of  deceleration. 
The  velocity  at  the  end  of  the 
deceleration. 
The  time  the  brake  is  not  in  use. 

Third  Card  (5110)  -  Thermal  analysis  information. 

a.  Cols  1-10   :  Number  of  nodes  (NEL)  . 

b.  Cols  11-20  :  An  integer  number  that  controls  the 
amount  of  printout  when  detailed 
output  is  required  during  the 
vehicle  deceleration.  The  amount  of 
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b. 

Cols  11-20 

c. 

Cols  21-30 

d. 

Cols  31-40 

lines  written,  depends  on  the 
stopping  time  and  time  increment. 
(NWR)  . 

c.  Cols  21-30  :  An  integar  number  that  controls 

the  amount  of  printout  when 
detailed  output  of  the  temperatures 
is  required  during  the  period  that 
the  vehicle  is  not  in  use.  The 
amount  of  lines  written  depends  on 
the  period  length  that  the  brakes 
are  not  used  (NWRA)  . 

d.  Cols  31-40  :  An  integer  number  that  controls  the 

amount  of  printout  when  detailed 
output  of  the  temperatures  are 
required  during  the  period  of 
constant  heat  generation.  (NWRQ) . 

e.  Cols  41-50  :  Number  of  braking  shoes  in  the 

machine. 
Fourth  Card  (7F10.0)  -  Brake  dimensions. 
a.   Cols  1-10   :  Inside  drum  radius  (RI) . 

Drum  width  (WDTH) . 

Drum  thickness  (THK) . 

Ratio  of  distance  from  pivot  to 

center  of  rotation  and  inside 

radius  (RD)  . 

e.  Cols  41-50  :  Drum  density  (RHO) . 

f.  Cols  51-60  :  Angle  between  hinged  pin  and  the 

begining  of  the  lining  (TETA1)  . 

g.  Cols  61-70  :  Angle  between  hinged  pin  and  the 

end  of  the  lining  (TETA2)  . 
Fifth  Card  (7F10.0)  -  Thermal  and  friction 

information, 
a.   Cols  1-10   :  Heat  conduction  coefficient  (K) . 

(real  number)  . 
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b. 

Cols 

11-20 

c. 

Cols 

21-30 

d. 

Cols 

31-40 

b. 

Cols 

11-20 

c. 

Cols 

21-30 

d. 

Cols 

31-40 

9. 

Cols 

41-50 

f. 

Cols 

51-60 

g- 

Cols 

61-70 

Heat  convection  coefficient  (HCNV) . 
Specific  heat  of  the  drum  (SPHT) . 
Max.  temperature  difference  between 
cold  friction  coefficient  and  hot 
friction  coefficient  (DFTM) . 
Cold  friction  coefficient  (CMIV) . 
Hot  friction  coefficient  (HSIV)  . 
Initial  temperature  (TINI) . 
Sixth  Card  (2F10.0)  -  Machine  information. 

a.  Cols  1-10   :  Vehicles  weight  (W) . 

b.  Cols  11-20  :  Wheel  radius  (RTIEE) . 
Seventh  Card  (5F10.0)  -  Analysis  constants. 

a.  Cols  1-10   :  Maximum  pressure  between  lining 

and  drum  (PRSA)  . 

b.  Cols  11-20  :  Constant  3.1415927 

c.  Cols  21-30  :  Gravitational  constant  (G) . 

d.  Cols  31-40  :  Increment  of  time  (ZHAN) . 

e.  Cols  41-50  :  The  angle  of  maximum  pressure 

(TETAA) . 
Eight  Card  (110.0)  -  Print  control. 

Cols  1-10   :  An  integer  number  can  be  zero  or  1. 
If  zero  (or  a  blank  card)  -  only 
the  final  results  are  printed. 
If  1-  the  temperature  at  time 
increments  are  printed. 
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Fig.  1   Dynamic  Representation  of  a  Brake  or  Clutch 
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'ig.  2   Brake  Assembly 
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Fig.  3   Finite  Difference  Model 
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